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Abstract We consider circular motion of a heavy object in an atomic Bose-Einstein 
condensate (BEC) at T = OK. Even if the linear velocity of the object is smaller 
than the Landau critical velocity, the object may radiate quasiparticles and thus 
experience the quantum friction. The radiation process is similar to Zel'dovich- 
Starobinskii (ZS) effect - the radiation by a rotating black hole. This analogy 
emerges when one introduces the effective acoustic metric for quasiparticles. In 
the rotating frame this metric has an ergosurface, which is similar to the ergosur- 
face in the metric of a rotating black hole. In a finite size BEC, the quasiparticle 
creation takes place when the ergosurface is within the condensate and occurs via 
quantum tunneling from the object into the ergoregion. The dependence of the 
radiation rate on the position of the ergosurface is investigated. 
PACS numbers: 03.75.Nt, 03.75.Kk 04.62.+V 04.70.Dy 



1 Introduction 

There are many challenges to simulate the phenomena on black holes by using 
condensed matter systems 1 . One of them is to study an analogy between gravity 
and superfluidity, in which a superfluid ground state (superfluid vacuum) serves 
as the analog of the vacuum of relativistic quantum fields, while the flow of the 
superfluid liquid imitates the metric field acting on 'relativistic' quasiparticles 
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(phonons, ripplons, fermionic excitations, etc.)^. The radiation from the black 
holes can be tested in this analogy, for example the Hawking radiation 1 . In this 
work we discuss how using atomic Bose-Einstein condensates one can simulate 
the Zel'dovich-Starobinskii (ZS) effect 3 '-. - the amplification and spontaneous 
emission of electromagnetic waves by a body rotating in the quantum vacuum, 
or by a rotating black hole (Kerr-BH). Radiation of electromagnetic mode occurs 
if its frequency ft) and its angular momentum quantum number / satisfy the condi- 
tion ft) IQ < 0, where Q is the angular velocity of the rotating body. This corre- 
sponds to the negative energy of the mode in the frame co-rotating with the body. 
In the Kerr black hole the negative energy states exist in the ergoregion, where the 
speed of the frame dragging exceeds that of light. In the ergoregion goo < 5 . The 
boundary of the ergoregion goo = 0, where the speed of frame dragging is equal 
to the speed of light, is called the ergosurface, 

2 The Zel'dovich-Starobinskii effect in superfluid vacuum 

Phonon excitation in a superfluid moving with velocity \ s has the Doppler shifted 
energy e(p) = cp + \ s ■ p. From this dispersion relation it follows that phonon 
trajectory in moving liquid obeys the null-geodesic equation ds = 0, where ds is 
the interval determined by the acoustic metric introduced by Unruh 6 

ds 2 = g liV dx tl dx v = c 2 dt 2 - (dr - \ s dt) 2 . (1) 

In superfluids the speed of light and the frame dragging velocity are played by 
sound velocity c and superfluid velocity \ s correspondingly. Then in the ergore- 
gion the superfluid velocity satisfies the condition v 2 > c 2 , while the ergosurface, 
on which goo = for the acoustic metric, is defined as a surface with v 2 = c 2 . 

Now let us consider an small object moving in a stationary superfluid along the 
circular trajectory with radius p = p t,j and with linear velocity v Q bj < c. Since in 
BEC the Landau critical velocity for creation of quasiparticles is VLandau = c > naive 
consideration suggests that the object is not radiating. However, radiation occurs 
due to Zel'dovich-Starobinskii effect. To see that let us choose the reference frame 
rotating with the angular velocity Q — v bj / p bj • In this frame the object is station- 
ary, while the superfluid moves with velocity field v s = — Qp§, where we use the 
cylindrical coordinate r = (p, <j),z). At p > c/Q, the superfluid velocity exceeds 
the speed of sound and goo < 0, which demonstrates that p erg = c/Q marks the 
position of the ergorsurface. In the ergoregion, i.e. at p > c/Q, the energy of some 
phonons becomes negative. These are phonons that satisfy condition ft) — IQ < 0: 

e(p) =cp + \ s p = cp-pp^Q — h(co-lQ) < , (2) 

It is energetically favorable to create these phonon modes in the ergoregion. How- 
ever, the source of the radiation - the moving body - is away from the ergoregion: 
Pobj = v bj /Q < Perg = c/Q.ln the process of radiation the phonon mode created 
by the source at p = p t,j is tunneling through the barrier of the positive energy 
region p e r g > p > p bj to tne ergoregion at p > p er g— • In the infinite system the 
radiation takes place for arbitrary Q, and thus for arbitrary velocity of the circulat- 
ing object, though for small v t>j c the ergoregion is far from the object and the 
probability of quantum tunneling is suppressed by factor (Qp^j/c) 21 = (v bj/c) 2 '. 
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3 The Zel'dovich-Starobinskii effect in atomic Bose-Einstein condensates 



Now let us consider how the ZS effect is modified in a finite system - in an atomic 
BEC trapped by an axisymmetrical potential at T = OK. The condensate is ini- 
tially at rest in the laboratory frame. The moving heavy body is reproduced by 
the additional potential which makes a small dip in the condensate at p = p t,j 
(see FigQ] left) and this dip is circulating with linear velocity v t>j < c. As in case 
of the moving body the circulating dip serves as an environment which provides 
us with the preferred reference frame. This is the frame rotating with frequency 
& — Vobj/Pobj m which the trapping potential and the dip are time-independent. 
The main difference from the infinite system discussed above is that in the finite 
system the negative energy states in the rotating frame appear only after Q reaches 
some critical value Q c . As a result the radiation occurs only after the velocity of 
the object reaches v c = £2 c p t,j. Thus in the finite system the critical velocity of 
the circulating object is non-zero, but it is much smaller than the Landau critical 
velocity c if the object is close to the center of the condensate. 
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Fig. 1 The plots of the density no(r) (a) and phase 0(r) (b) of the condensate with the object 
circulating with the angular velocity Q > Cl c . (a) The object potential makes a density dip. (b) 
The phase is almost uniform: its variation is on the order of 10~ 2 . Here the length x and y are 
scaled with b = ^Jh/mOi, and the density n$ of BEC is normalized: / n^dr = 1. 



The states with Q < Q c are frictionless: they correspond to the local min- 
ima of the energy E — L z £2, where L z is the angular momentum of the superfluid 
(Fig|2}(a)). In these states, all the quasiparticles have positive energy in the ro- 
tating frame; h((0 — IQ) > 0. At Q > Q c the stability is lost, the local minimum 
transforms to the saddle point (Figf2}(b)), and the moving object would radiate 
quasiparticles with negative energy h((0 — IQ) < 0. 

Let us illustrate this using numerical simulations of the BEC. In the absence of 
quasiparticles, BEC at T = can be described by the macroscopic wave function 
<£(!■) = \J no (r) exp [id (r)] , which satisfies the Gross-Pitaevskii (GP) equation 

V&=^-^V 2 + V tmp + V ohj +g\<P\ 2 -nL)jcp . (3) 

Here L z = —ih(r x V) z ; g = Aittra/m; in and a are the atomic mass and s-wave 
scattering length, V tlap and V (,j refer to the trapping and the object potential, re- 
spectively. We shall consider the two-dimensional system with a harmonic poten- 
tial Vti a p = 1 /2m(6{ -p , and the object potential Vobj is assumed to be localized and 
repulsive. The frictionless states below Q c are obtained using the imaginary time 
propagation (ITP) of the Gross-Pitaevskii energy. The states above Q c are also 
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(a) E{Q<Q C ) 



(b) 
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Fig. 2 The schematic diagram of the energy of the BEC in the frame co-rotating with the circu- 
lating body, (a) The state at Q < Q c corresponds to the local minimum and keeps the frictionless 
circulation of the object. The energy gap which stabilizes the frictionless state is reduced as Q. is 
increased, (b) The state at Q. > Q c is energetically unstable and is destroyed by the spontaneous 
radiation of quasiparticles with negative energies in the rotating frame. 



obtained with the ITP if the time of development of instability is long enough. In 
FigQjthe state is shown with Q. = 0.56 x a\ > £2 C , g = 200 x h(Ob 2 /N2D, where 
N2D ls the particle number per unit length along the z axis and b = \Jfijm(i\. 
The object potential has a Gaussian form V bj = Vb ex p[~ ( r — r obj) 2 /w 2 ], where 
V = 0.05 x mtfR 2 T , r obj = (-0.1 x R T ,0), w = 0.05 x R T and R T = 3.99 x b is 
the Thomas-Fermi radius of the condensate. 

Since the object potential is weak, the circulating object makes a small density 
dip at r = r Q bj (Figfl]-(a)), which practically does not disturb the superfluid veloc- 
ity (7i/m)V0(r): it is almost zero everywhere in the laboratory frame (FigQ]-(b)). 
The role of the dip is to provide the perturbation which violates the cylindrical 
symmetry and induces the nonzero matrix element for the transition between the 
state without quasiparticles and the states with quasiparticles; the dip also pro- 
vides the distinguished reference frame which in the absence of other interactions 
serves as the frame of the environment. 



4 Critical velocity of circulating body 

Let us now find the critical velocity Q c , at which the negative energy phonon 
mode appears in the condensate and thus the radiation becomes energetically pos- 
sible. The spectrum of quasiparticles are obtained by coupling the GP eq.© and 
the Bogoliubov de Gennes (BdG) equation. For the calculation of the spectrum 
we may neglect the small deformation of the condensate caused by the object 
potential and consider the axisymmetric BEC ^(r) = ^/no(p). Then the Bose 
operator is written as 'f'(r) = <P(p) + Y!ki[ u ki&ki ~ v u*&li}^ where u k i(r,t) = 
a^(p)exp[/(/(/) — £ k it /%)] and Vki(r,t) = Vki(p)exp[i(l^) — £ u t /%)]. The numbers 
k and / refer to the radial and angular quantum number of quasiparticles and the 
sum Y! ls taken without a trivial zero-energy mode (k, I) = (0,0). As a result the 
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BdG eq. is reduced to 8 
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-h 2 /2m[(d/dp) 2 + (l/p)(d/dp) - l 2 /p 2 ] + y tap +2g\<P\ 2 - 
■ g<P 2 . The normalized condition of and v\\ is / t/r ( | Mjy | 2 — | v# 
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1. It can be shown that the values u# and Vjy are independent of Q. but the energy 
spectrum depends as £/t/(I2) = £,«(0) — hlQ. 
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Fig. 3 (a) Critical velocity Cl^ for nucleation of modes (k = 0,1) and (k = 1,1) as a function 
of the orbital quantum number /. The critical velocity Sl c is given by the lowest of them: £2 C = 
CI-] ~ 0.508tt)t. (b) Position of the ergosurface p elg in a semiclassical description as a function of 
12. (c) Probability density «o; of the modes (0,/); «oo = «o refers to the condensate density, (d) 
The radiation rate of the mode / by the body is proportional to w'o/(p bj) = «o(pobj)' ! o;(Pobj) at 
the position of the body, wqi is plotted as a function of the position of the body. 



The mode (k,l) can be radiated if the rotation velocity Q reaches the value 
Qki, which satisfies condition e^;(i2<.;) = 0. This value Q.^ is plotted in Figj3]-(a) 
for g — 200 x h(0tb 2 /N2D- Modes with (0, 1) are most important, since modes with 
higher k have higher Q^. From Fig.[3]-(a) it follows that the mode (k,l) = (0,7) 
has the lowest value of the critical velocity; it thus determines the global critical 
velocity at which radiation becomes possible: Q c = Qqj ~ 0.508<i)f 



5 Radiation by circulating body 

Since quantum numbers of relevant modes are small, the classical notion of the er- 
goregion is not well defined. Nevertheless, ergosurface may be crudely introduced 
as the surface p erg = c(p erg ) where c(p) = \J grio(p) / in is the local speed of 
sound (Figj3]-(b)). The position of the ergosurface for Q = Q c is p el -g(£2 c ) — 3.2, 
i.e. at the periphery of the condensate and far from the circulating object. This 
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confirms the general picture that an object circulating with the subsonic velocity 
spontaneously radiates sound modes which are formed far from the object. 

This is also supported by exact quantum mechanical consideration. If the cir- 
culating object is small the matrix element of transition is concentrated locally 
at the position of the object r = r (,j, and the radiation rate of the mode (k,l) is 
proportional to the condensate density «o(robj) at the position of the object and to 
the probability density «H( r obj) of the mode also at the position of the object, and 
thus it is determined by the product w^Pobj) = «o(Pobj W/(Pobj)- 

The probability density of the (0,/) mode «o/(p) = l M o/(p)| 2 + l v o/(p)| 2 is 
shown in FigOKc). The 'penetration' of the mode into the central region of the 
condensate is monotonically decreasing as p is reduced. In terms of the semiclassi- 
cal description, this reduction can be interpreted as the tunneling of quasiparticles 
into the energetically forbidden region 2 . The quantity Wki(p bj) = n o(Pobj) H A-z(Pobj) 
which determines the radiation rate is shown in Figf3}(d). It demonstrates that the 
radiation is enhanced as the object becomes closer to the ergosurface. 



6 Conclusions 

Spontaneous radiation is possible by an external object circulating in BEC with 
linear velocity much smaller than the Landau critical velocity. The origin of radi- 
ation is similar to the ZS radiation from the Kerr black hole. It can be described 
in terms of quantum tunneling through the energetically forbidden region between 
the object and the ergoregion. The radiation rate depends on rotation speed and on 
the distance between the object and the ergosurface. The role of the object might 
be played by the moving optical potential made by the blue or red detuned laser. 
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